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Abstract 

The modeling of the elastic properties of granular or nanoscale systems requires 
the foundations of the theory of elasticity to be revisited, as one explores scales at 
which this theory may no longer hold. The only cases for which a microscopic jus- 
tification of elasticity exists are (nearly) uniformly strained lattices. A microscopic 
theory of elasticity, as well as simulations, reveal that standard continuum elasticity 
applies only at sufficiently large scales (typically 100 particle diameters). Interest- 
ingly, force chains, which have been observed in experiments on granular systems, 
and attributed to non-elastic effects, are shown to exist in systems composed of 
harmonically interacting constituents. The corresponding stress field, which is a 
continuum mechanical (averaged) entity, exhibits no chain structures even at near- 
microscopic resolutions, but it does reflect macroscopic anisotropy, when present. 
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1 Introduction 



It is quite surprising that the existing microscopic justification of the time- 
honored theory of elasticity, which has been thoroughly researched in a variety 
of disciplines, is limited to lattice atomic configurations [1]. Classical contin- 
uum elasticity theory has been applied to a large variety of systems, including 
granular materials [2] . In recent years the same theory has been applied for the 
description of elastic properties of micro- and nano-scale systems (e.g., [3]). It 
is a-priori unclear whether this theory applies at such small scales. 
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The study presented below shows that the jTistification of elastic theory based 
on a 'microscopic' picture is not entirely straightforward. In fact, we show 
that continuum elasticity is limited to sufficiently large scales. One of the other 
results presented in this paper is that 'force chains' in granular systems, which 
have been attributed by some authors to non-elastic effects [4] , are reproduced 
in simulations of elastic systems, in good agreement with experiments. 



2 Coarse Graining and Constitutive Relations 

2. 1 Preliminaries 

Below the term 'particles' is taken to mean grains in a granular system or 
atoms in a solid. Classical mechanics is assumed throughout this paper. Fol- 
lowing [5] define, for a system of particles, the coarse grained mass and mo- 
mentum densities, at the position r and time by p(r, t) = I]jmj0[r — rj(t)] 
and p(r,t) = I]j TTijVj(t)0[r — rj(t)], respectively, where {rj(t); Vj(t); mj} are 
the positions, velocities and masses of the particles, indexed by i, and 0(R) is 
a normalized non-negative coarse graining function (with a single maximum 
at R = 0) of width A, the coarse graining scale. The velocity field is defined 
by V(r, = p(r, f)/p(r, i). Upon taking the time derivatives of these coarse 
grained fields and performing straightforward algebraic manipulations [5] one 
obtains the equation of continuity, as well as the momentum conservation 
equation: p„(r, t) = -g^ [p(r, t)Va(r, t)V^(r, t) - (Ta^(r, t)], where Greek in- 
dices denote Cartesian coordinates. The stress tensor, aap, is composed of 
a kinetic contribution, irrelevant for our purposes (we consider quasi-static 
deformations), and a "contact" contribution. Neglecting the former, one ob- 
tains: 

1 } 
o-a/3(r,t) = -- fija{t)rij^{t) j ds4)[r - ri{t) + srij{t)]. (1) 

where fijait) is the a-th component of the force exerted on particle i by particle 
j {j ^ i) at time t (assuming pairwise interactions) and 

2.2 Displacement and Strain 

Following elementary continuum mechanics, the velocity of a material particle 
whose initial (Lagrangian) coordinate is R satisfies: V^'^(R, t) ~ du^^ (R, t)/dt, 
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where u'"'^(R, t) = r(R, t) — R is the (Lagrangian) displacement field. There- 
fore: u^''(R, t) = J^Y^^{'R,t')dt'. Using the definitions presented in Sec. 2.1, 
it follows that: 

i Sim>[r(R,«')-r^(«')l ^ ' 



Noting that lij = v^, where is the displacement of particle i, and invoking 
integration by parts in Eq. (2), one obtains: 



where the result has been translated to the Eulerian representation. The cor- 
rection terms can be shown to beget non-linear terms in the constitutive rela- 
tions. The linear strain is given by eaf3{r,t) = | ^"g^^^'*^ -|- . Note that 
the exact displacement and strain obtained here are different from the com- 
monly defined mean field strain, cf. e.g. [6]. Mean field theories are based 
on the assumption that the relative particle displacements are described by 
the macroscopic strain, Uija{r, t) — ea/sRijis, where Uij = Uj — u^. A suggested 
improvement is the "best fit" strain [7]. These mean field approaches are in- 
accurate and often lead to incorrect constitutive relations [8]. 



2.3 Stress-Strain Relation 



The microscopic expressions for the stress [Eq. (1)] and the strain are not 
manifestly proportional. Therefore elasticity is not a-priori obvious. To see 
how it still comes about, consider the case of harmonic (and local) pairwise 
particle potentials Uij = \Kij (|rjj| — /y)^. The forces are then given, to linear 

order in 5vij = Uij, by fij = —Kij (r'fj ■ 6rij^ rfj. The stress [Eq. (1)] is then 
given by a^;^(r, t) = ^ J2ij Kijr'^j^Srij^f'^j^r^jf^(j){r - r^), up to nonhnear terms 
in {6rij}. The superscript denotes the reference configuration. 

Consider a volume Q which is much larger than the coarse graining scale, 
A, and let r be an interior point of Q which is 'far' from its boundary. Let 
upper case Latin indices denote the particles in the exterior of ^2 which in- 
teract with particles inside il. Since the considered system is linear, there 
exists a Green's function G such that — Giaj/3Ujp for i e 0. Hence, 
Uija = {GiaJiS — Gjajis) Uj3- Whcu all Uj arc equal (rigid translation), Ujj = 0. 
Hence, Uija = {Giaj/s — Gj^jp) [mj/j — Uisi^)]- Define Lijajp = Giajp — Gjajfs- 
One then obtains Uija = Lijajp [upirj) - ^/^(r)] + Lijajp [ujp - «/3(rj)], where 
ujp — upir j) is Sk fluctuating displacement. The sum over J in the second term 
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can be shown to be subdominant when the hnear size of the volume Q greatly 
exceeds the coarse graining scale. Furthermore, to lowest order in a gradient 
expansion, M^(rj) — M^(r) ~ {rj^ — r^). Substituting this result in cr^^^, 

and employing rotational symmetry, one obtains: 

(^a/sir) ^ ^e^^{r) ^ Ki^L^^^j^ (r°^ - r^) r° ^0(r - r°). (4) 

Thus linear elasticity is valid when |AVq,V/3u| -C 1 and ||e|| -C 1. Note that the 
elastic moduli depend, in principle, on the position as well as the resolution 
(through the coarse graining function 0). 



3 Numerical Results 

In order to study the crossover from microelasticity to macroelasticity [9], 
we have considered a two dimensional (2D) rectangular slab, composed of 
uniform disks (of diameter d) whose centers are positioned on a triangular lat- 
tice, with nearest neighbors coupled by uniform linear springs (of rest length 
d), and a three dimensional (3D) slab, composed of uniform spheres (of di- 
ameter d) with springs coupling nearest-neighbors (spring constant Ki) and 
next-nearest-neighbors (^^2)- Both systems can be shown to correspond to 
isotropic elastic media in the continuum limit [1] (for the 3D system, only if 
Ki — K2). The particles at the bottom layer are attached to a rigid floor, 
and a downward force is applied to the central particle in the top layer. This 
is intended for comparison with 2D [10] and 3D [11] granular experiments 
(even though these models do not directly describe the interaction of gran- 
ular particles and are primarily intended for studying the crossover between 
microelasticity and macroelasticity, they closely reproduce some of the ex- 
perimental results). Fig. 1 presents a comparison between the vertical stress 
exerted on the floor of the system for different slab heights (numbers of layers 
of particles) and continuum elastic solutions, for both systems {Ki = K2 for 
the 3D system). The convergence to the appropriate (rough rigid support) 
elastic stress distribution for a sufficient number of layers is evident. 

The observation of "force chains" in experiments [10] and simulations [12] has 
been interpreted as evidence for non-elastic effects [4]. However, in a plot of 
the interparticle forces in the 2D system described above (Fig. 2, left) one 
can identify "force chains" though the microscopic interactions are strictly 
linear, and the large scale behavior corresponds to isotropic elasticity. The 
corresponding distribution of interparticle forces [9] closely resembles the ex- 
perimental findings of Geng et al. [10]. The zz component of the stress field, 
calculated using Eq. (1) with a Gaussian coarse graining function [5] of width d 
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Fig. 1. The vertical stress at the bottom of a 2D triangular lattice (left) and a 3D 
cubic lattice (right), compared to continuum elastic solutions. The results are scaled 
by the slab height h. 

(i.e., with fine resolution), is shown in Fig. 2, on the right. The force chains are 
not evident any more. The force chains result from the microscopic anisotropy 
dictated by the fact that a single particle is never subject to an isotropic force 
distribution. However, even at small, but finite, spatial resolution, the calcu- 
lation of the stress tensor inherently involves averaging over the forces acting 
on the particles, so that this anisotropy does not appear in the stress field. 
Similar results [9] are obtained for random 2D systems. 




Fig. 2. Left: Force chains in a 2D triangular lattice. A vertical force is applied to the 
center particle in the top layer. The line widths are proportional to the forces. Only 
the central region of the lattice is shown (the lattice comprises 15 x 41 particles). 
Right: Contour plot of the vertical stress component, azz, in the same region. 

A more realistic model with unilateral ( "one-sided" ) springs, active only upon 
compression, yields results which are even closer to the experiments [9] . In this 
case, some of the springs are disconnected (as obtained in [13] for a pile geom- 
etry) when the external force is applied. In particular, the horizontal springs 
in a triangular region below the point of application of the force are severed, 
and the corresponding macroscopic properties are anisotropic. The resulting 
stress distribution [9] is similar to that predicted by hyperbolic models of force 
propagation [4], but it is a result of an elastic model. 
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